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ABSTRACT 



Context. The features of pulsed y-ray emission from classical and millisecond pulsars indicate that the high energy radiation processes 

in their inner magnetospheres occur in a similar way. In the last decade several TeV y-ray nebulae have been discovered around 

classical pulsars. The above facts suggest that y-rays should be produced also in the surroundings of millisecond pulsars. 

Aims. We discuss a model for the bow shock nebula around the well known Black Widow binary system containing the millisecond 

pulsar B 1957+20. This model predicts the existence of a synchrotron X-ray and inverse Compton y-ray nebula around this system. 

We want to find out whether y-ray emission from the nebula around B 1957+20 could be detected by the future and present Cherenkov 

telescopes. 

Methods. Using the Monte Carlo method we followed the propagation of relativistic electrons in the vicinity of the pulsar. We 
calculated the very high energy radiation produced by them in the synchrotron process and the inverse Compton scattering of the 
Microwave Background Radiation and of the infrared radiation from the galactic disk. We also computed the X-ray emission produced 
by the electrons in the synchrotron process. 

Results. We show that the hard X-ray tail emission observed from the vicinity of B 1957+20 can be explained by our model. Moreover, 
we predict that the TeV y-ray emission produced by the electrons in the inverse Compton process should be detectable by the future 
Cherenkov Telescope Array and possibly by the long term observations with the present Cherenkov arrays such as MAGIC and 
VERITAS. The y-ray emission from B 1957+20 is expected to be extended, inhomogeneous, and shifted from the present location of 
the binary system by a distance comparable to the radius of the nebula. 

Key words, pulsars: general — stars: binaries: close — radiation mechanisms: non-thermal — gamma-rays: general 



> 



X 



1. Introduction 

PSR B 1957+20 was the first millisecond pulsar (MSP) discov- 
ered within the binary system belonging to the class of Black 
Widows (Fruchter et al. 1988). This pulsar has a very small mass 
companion (~ 0.022Mq, van Paradijs et al. 1988) which evapo- 
rates under the irradiation from the pulsar magnetosphere. The 
pulsar has the period of 1 .607 ms, the surface magnetic field of 
~ 10*^ G, and the rotational energy loss rate of 7 .5 x 10^"* erg s"' . 
The distance to the binary system is estimated on 2.5 kpc (from 
the model for Galactic electron density) consistent with the re- 
cently established lower limit ~2 kpc (van Kerkwijk et al. 201 1). 
The binary system is compact with the orbital radius of 1 .5 x 10" 
cm. The companion star has the radius ~ 10"^ cm and the surface 
temperature which varies between 2900 K for the unilluminated 
side to 8300 K for the illuminated side (Fruchter et al. 1995, 
Reynolds et al. 2007). Therefore, stellar radiation is not expected 
to create a very strong target for relativistic particles within the 
binary system. At present, the companion star loses mass at a 
rather low rate estimated on 10""^ M© yr"' (Takata et al. (2012). 

The importance of the high energy processes in the vicinity 
of PSR B 1957+20 has become clear with the discovery of an Ha 
emission nebula (KuUcami & Hester 1988). This emission is ex- 
pected to be produced in shocks formed in the interaction of the 
pulsar wind with the interstellar medium. A clear bow shock has 
been detected which apex is located at the distance of ~ 0.02 pc 
from the pulsar. The bow shock appears due to the motion of 



the binary system with the velocity 220 km s"' through the in- 
terstellar medium (Arzoumanian et al. 1994). The X-ray emis- 
sion has been also reported from the direction of the binary sys- 
tem in the observations of Chandra (Stappers et al. 2003, Huang 
& Becker 2007, Huang et al. 2012). This emission comes from 
the interior of the bow shock creating a tail behind the mov- 
ing binary system. The length of the tail is ~ 10'^ cm (Huang 
et al. 2012). The X-ray emission is well described by a single 
power law spectrum with the index in the range 2.3 - 2.6 de- 
pending on the absorption model. The extended X-ray feature 
has been interpreted as emission from energetic electrons which 
radiate on the crossing time scale of this region by the pulsar 
moving with velocity of 220 km s"' (Cheng et al. 2006). 

The Black Widow binary system containing B 1957 +20 was 
claimed in the past to be a GeV-TeV y-ray source (Brink et 
al. 1990). But this early report was not confirmed in the anal- 
ysis of the EGRET data (Buccheri et al. 1996). In fact, such high 
energy emission has been suspected akeady since the discovery 
of Black Widow pulsars as a result of either the acceleration of 
particles within the binary system or in the the shock waves of 
the pulsar wind (e.g. Arons & Tavani 1993, Cheng et al. 2006, 
Takata et al. 2012). Recently, a pulsed GeV emission from the 
pulsar B 1957+20 has been discovered by Fermi (Guillemot et 
al. 2012). The pulsed spectrum is flat above 0. 1 GeV (spectral in- 
dex close to 2) and extends up to ~ 4 GeV. The phasogram (light 
curve folded with the period of the pulsar) shows two well sepa- 
rated peaks. Such structure is also common in the case of classi- 
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Fig. 1. Schematic representation of the bow shock nebula around 
binary system containing the millisecond pulsar B 1957+20. The 
bow shock is created due to the motion of the binary system 
through the interstellar space with the velocity of ~ 220 km s~' . 
Relativistic electrons with the Lorentz factors ye are accelerated 
by the pulsar itself or by the shocks due to the pulsar wind inter- 
actions. The electrons are coUimated by the bow shock in the di- 
rection opposite to the motion of the binary system. These elec- 
trons comptonize the Microwave Background Radiation (MBR) 
and the infrared radiation (INF) from the galactic disk. As a re- 
sult y-ray photons are produced (tagged as Ey) at the region be- 
hind the pulsar. 



cal pulsars. Therefore, it is expected that the radiation processes 
in the inner magnetosphere of the millisecond pulsar B 1957 +20 
are similar to those occurring in the case of classical pulsars. 
This strongly indicate that also processes of acceleration of par- 
ticles in the pulsar wind are expected to occur similarly. Very 
recently, Wu et al. (2012) reports detection of orbital modula- 
tion of the y-ray emission at energies above ~2.7 GeV from the 
Black Widow pulsar PSR B 1957+20. This emission is expected 
to be produced by electrons from the pulsar wind which comp- 
tonize stellar radiation. 

We investigate the radiation processes in the supposed pul- 
sar wind nebula around the binary system containing PSR 
B 1957+20. The synchrotron X-ray and inverse Compton (IC) 
y-ray emission is calculated from such nebula for the range of 
likely parameters. Based on the comparison of the calculated 
synchrotron spectrum with the observed X-ray emission we con- 
clude on the detectability of the TeV y-ray emission from the 
bow shock nebula surrounding PSR B 1957+20. 



2. The nebula around binary system downstream of 
the bow shock 

Since the proprieties of high energy y-ray emission from the 
milUsecond pulsars and classical radio pulsars are surprisingly 
similar (see the first pulsar catalogue, Abdo et al. 2010), it 
seems clear that the processes occurring in their inner magne- 
tospheres are these same. Therefore, millisecond pulsars should 



also produce relativistic pulsar winds with the parameters simi- 
lar to those observed around classical pulsars. However, nebulae 
around MSPs are expected to have a very complicated structure 
(and also other proprieties) since many MSPs form compact bi- 
nary systems which additionally move in the interstellar space 
with large velocities. In fact, this is the case of the binary system 
PSR B 1957+20. The pulsar wind around B1957+20 is expected 
to interact with the induced wind of the low mass companion 
star within a small solid angle, of the order of -0.01 sr, corre- 
sponding to eclipse time of the pulsar radio emission by the wind 
of the companion star (Fruchter et al. 1988). Therefore, most of 
the pulsar wind is expected to escape unaffected from the binary 
system. Due to the fast velocity of the binary system, the pulsar 
wind has to interact with the interstellar medium creating a bow 
shock. Such bow shock has been detected in Ha emission in the 
case of B 1957+20. The distance of the apex of the bow shock to 
the pulsar is estimated on ~ 0.02 pc (Kulkami & Hester 1988). 
This bow shock confines the pulsar wind at least in the direc- 
tion of pulsar's motion. Relativistic electrons in the wind can 
diffuse mainly in the direction opposite to the pulsar's motion. 
Cheng et al. (2006) have proposed that the synchrotron radia- 
tion from such ultrarelativistic electrons is responsible for the 
observed X-ray tail extending along the axis of the bow shock 
in the direction opposite to the pulsar velocity. We intend to per- 
form calculations of the Inverse Compton (IC) y-ray emission 
from such relativistic electrons applying general scenario pro- 
posed by Cheng et al. (2006). These authors argue that efficient 
synchrotron emission by electrons can occur on the dynamical 
time scale of the pulsar crossing the length of the tail estimated 
on ~10^^ cm (Huang et al. 2012). This dynamical time scale is 
equal to, 

Tdyn=/?/Vbm~ 1.5X10"/?1 S, (1) 

where the velocity of the binary system is Vbin = 220 km s"^ and 
the length of the tail is = IR\ pc. Applying the observed length 
of the synchrotron emission, we can estimate the optical depth 
for electrons on the IC scattering of the Microwave Background 
Radiation (MBR) in the Thomson regime (true for electrons with 
energies below ~ 100 TeV) on, r = CTdyniMBRO'T ~ 0.35, 
where c is the velocity of light, crj is the Thomson cross section, 
and hmbr is the photon density of the Microwave Background 
Radiation. Note however that electrons cool only partially in the 
region of observed X-ray emission. Many of them escape from 
this region but continue to interact with the MBR and other soft 
photon field, producing high energy y-rays. Therefore, we ex- 
pect the appearance of the y-ray nebula in the vicinity of the 
Black Widow binary pulsar. This nebula should be shifted in re- 
spect to the observed location of the binary system in the direc- 
tion opposite to the pulsar's motion. 

On the other hand, the energy loss time scale of electrons on 
the IC scattering in the Thomson regime is, 



OTeC ye/(4cf/radO'Tye/3) S, 



(2) 



'IC 

where Umd is the energy density of the soft radiation field equal 
to 0.3 eV cm"^ for the Microwave Background Radiation (MBR) 
and to ~ 1.5 eV cm"^ for the infrared radiation with characteris- 
tic energies ~ 0.01 eV, produced in the galactic disk (e.g. see the 
values calculated in Hui et al. (2011) based on the GALPROP 
code developed by Strong & Moskalenko 1998), and me is the 
rest mass of an electron. For these energy densities we obtain the 
energy loss time scales of the order of ~ 6.3 x 10^^/ye s for 
the MBR and ~ 1.3 x lO'^/ye s for the infrared radiation, where 
ye is the Lorentz factor of the electrons. In order to cool the elec- 
trons efficiently on the IC process during the dynamical time of 
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the moving pulsar, the emission region should have the diameter 
of the order of R ^ 8.5 x 10^/ye pc. For example, in the region 
of 10 pc, electrons with energies larger than ~ 4 TeV (but below 
~ 100 TeV since the electrons have to interact in the Thomson 
regime) should be able to produce efficiently y-rays in the IC 
process by scattering infrared photons from the galactic disk. 
Note that, the region of the y-ray production in the IC process 
should be clearly shifted from the pulsar position in the direc- 
tion of the observed tail X-ray emission. This region should be 
also inhomogeneous with higher energy y-rays produced closer 
to the pulsar. 

In the above estimates we neglected the energy density of 
stellar photons, in respect to the MBR and infrared radiation 
at the region of the acceleration of electrons (the shock in the 
pulsar wind). In fact, the energy density of stellar photons de- 
pends on the distance from the star as f/* « 4.5 x 10~^/Z)jg eV 
cm"^, where the distance from the star is D = 10'**Z)i8 cm. It 
is assumed that the companion star in the binary system PSR 
1957-1-20 has the radius 10'" cm and most of its surface has tem- 
perature close to ~3000 K (Fruchter et al. 1995). For these pa- 
rameters, the electron energy losses are dominated by scattering 
of the infrared photons for distances above ~ 5 x lO'"" cm. Note 
also that the scattering of the optical photons from the star occurs 
in the Klein-Nishina regime for electrons with energies above 
~ 100 GeV. Therefore, the effects of scattering stellar radiation 
by the TeV electrons can be safely neglected. 

The region of the y-ray production can be also affected by 
the diffusion of the electrons in the pulsar wind downstream of 
the pulsar wind shock. We estimate the diffusion distance of the 
electrons, as a function of their energy, and compare it with the 
time scale corresponding to the dynamical motion of the pul- 
sar. For the Bohm diffusion approximation, the diffusion dis- 
tance is Rdif = yllD^iit, where Ddif = c/?l/3 is the diffusion 
coefficient, /?l is the Larmor radius of electrons, B is the mag- 
netic field strength in the considered region, and t is the diffusion 
time. If B is fixed on 1/iG, then Ddif » 1.5 x lO'^ye cm^s ' and 
/?dif ~ 5.5 X 10^ VXe^ cm. The spread of the emission region due 
to the diffusion process is smaller than that one due to the motion 
of the pulsar, i.e. Rm < /?dyn - Vpuif, for the following condition 
t > 6.2 X lO^'ye s. We compare this condition with the energy 
loss time scale on the IC process in the Thomson regime (see 
Eq. 2 and estimates below). It is found that electrons with ener- 
gies below ~7 TeV lose energy on production of y-rays when the 
ballistic motion of the binary system determines the morphology 
of the y-ray source. We conclude that depending on the electron 
energy, the dimension of the y-ray source is determined either 
by the motion of the Black Widow binary system through the 
interstellar medium or by the diffusion process of the electrons. 



3. Relativistic electrons in nebula 

We estimate the magnetic field strength around the pulsar, above 
its hght cyhnder radius, by extrapolating it from the pulsar sur- 
face. The magnetic field strength is then given by. 



BiR) ^ 4.4 X 10-V/2fi8/(/'ms/?i8) G, 



(3) 



where R = 10'^/?i8 cm is the distance from the pulsar, Bns - 
10^ B^ G is the magnetic field strength on the neutron star sur- 
face, P = 10"^Pms s is the period of the millisecond pulsar, 
and o" is the magnetization parameter of the pulsar wind, cr has 
been estimated in the case of the Crab Nebula on 0.003 (de 
Jager & Harding 1992) and in the case of the Vela Nebula on 
~0.1 (Sefako & de Jager 2003). cr is expected to be in the range 



0.001 - 0.01 in the modeling of the Crab Nebula presented by 
Kennel & Coroniti (1984). The magnetic field given by Eq. 3, 
is expected to be enhanced at the shock region in the pulsar 
wind by a factor of ~ 3. Downstream of the shock, electrons 
are isotropized and start to radiate efficiently synchrotron radia- 
tion. Therefore, the magnetic field in the region downstream of 
the shock is an important factor which determines the diffusion 
of the relativistic electrons and production of the synchrotron ra- 
diation. The maximum energies to which the electrons can be 
accelerated in the pulsar shock region can be estimated from. 



cR,hB(R,h) 10Vi/2B,/p2 GeV. 



(4) 



Note that this simple formula gives the energies of electrons 
present in the Crab Nebula comparable to those expected from 
the modelhng of its multi-TeV y-ray spectrum (e.g. de Jager & 
Harding 1992). 

As noted above, Chandra has detected the tail behind the 
pulsar B1957H-20 in the energy range 0.3-8 keV (Huang et 
al. 2012). If this emission is due to the synchrotron process from 
the relativistic electrons, then the Lorentz factors of the electrons 
can be estimated from. 



e = m^c\BIB„)yl, 



(5) 



where e = 8 keV is the energy of synchrotron photons, B and 
Ba - 4.4 X lO'^ G are the magnetic field in the emission re- 
gion and the critical magnetic field strength. The inspection of 
the above equations allows us to conclude that the production 
of the synchrotron photons with ~10 keV energies is possi- 
ble provided that the Lorentz factors of electrons are at least 
ye = 2.2 X \Q^P^,R^^^I{cr^l'^B^^\ obtained by substitution of 
Eq. 3 into Eq. (5). Electrons are accelerated to such energies 
provided that the magnetic field is strong enough, i.e. the shock 
in the pulsar wind appears close to the pulsar. For the parame- 
ters of PSR B 1957-1-20, the distance of the shock has to be below 
/?i8 ^ 3.4 X lOV^''^, which equals to R^u ^ W - 10^" cm for 
a in the range 0.001 - 0.1. This condition is consistent with the 
observations of the PWNe around classical pulsars. For exam- 
ple, in the case of the Crab Nebula the shock is located at the 
distance of ~ 3 x 10" cm (Kennel & Coroniti 1984). 

It is not clear at present in what process electrons reach such 
large energies. This might be reconnection of the magnetic field 
or the shock acceleration mechanism. In the second case, the 
limit on the maximum energies of the electrons have to be con- 
sistent with the hmit due to the presence of the synchrotron en- 
ergy losses already during the acceleration process. The maxi- 
mum energies of the electrons, due to the saturation by the syn- 
chrotron energy losses, can be derived from the comparison of 
the electron acceleration time scale. 



Tacc ~ ^EJ(X-lB) S, 

with the synchrotron energy loss time scale, 



' syn 



Ee/^syn ~ 370/(B2£) s, 



(6) 



(7) 



where £syn = i4/3)cUBO-TEi/mi » 0.0027 B^E^ TeV/s, the ac- 
celeration efficiency is parametrised by the factor = W^x-i' 
and Ee is the electron energy in TeV. Energies of the electrons 
can not be larger than, 

P R^'^ 

2 X 10\-^)'/2 GeV^5xlO'^ GeV, (8) 

^ B o-^I^b]^^ 

o 

For the pulsar with the parameters of PSR B 1957+20, Zimax is 
lower than is^ax for the location of the shock at R\% > 0.05cr^'^, 
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which corresponds to /?sh > 1 .6x 10'^ cm for cr = 0. 1 . Therefore, 
we conclude that for the expected localizations of the shock in 
the nebula around the pulsar B 1957 +20 (above ~ 10'^ cm), the 
synchrotron energy losses can not limit the acceleration process 
of the electrons below the maximum possible energies given by 
Eq.a 

4. Production of high energy radiation 

We calculate the y-ray spectra produced by relativistic elec- 
trons in the IC scattering of the MBR and the infrared radia- 
tion from the galactic disk. These electrons also produce syn- 
chrotron emission which can extend up X-ray energy range. It is 
commonly expected that electrons accelerated at the pulsar wind 
shock obtain the power law spectrum. We assume that this spec- 
trum has a lower energy cut-off at energies corresponding to the 
Lorentz factor of the pulsar wind, i.e is equal to a few times 
10^. In our calculations we fix this value on 3 TeV, in agreement 
with the modelling of the PWNe (Kennel & Coroniti 1984) and 
recent calculations of the spectra of the electrons leaving the in- 
ner magnetospheres of the millisecond pulsars in the frame of 
the pair starved polar cap model (e.g. Zajczyk et al. 2010). The 
electrons take a significant part of the energy lost by the mil- 
lisecond pulsar, which is of the order of ~ 10%. The spectrum 
of the electrons extends up to the maximum energy described in 
Sect. 3. These electrons are accelerated close to the pulsar wind 
shock and diffuse to the outer region creating a tail trailing be- 
hind the pulsar In this calculations we take the energy density 
of the infrared galactic disk emission equal to 1 .5 eV cm"^. It is 
assumed that the magnetic field is enhanced by a factor of 3 in 
the pulsar wind shock and at larger distances continue to drop 
according to Eq. 3 up to the minimum value Bmin. This mini- 
mum magnetic field strength can be even below the magnetic 
field strength in the interstellar space (of the order of ~2-6//G), 
since the volume of the pulsar wind nebula is not penetrated by 
the interstellar medium. 

We assume that electrons are injected at the distance of the 
shock from the pulsar, R^h- They slowly diffuse outward ac- 
cording to the Bohm diffusion model in a decreasing magnetic 
field. During the diffusion process, the electrons interact with 
the background radiation producing GeV-TeV y-rays in the IC 
process. We apply the Monte Carlo method in order to deter- 
mine the energy of the y-ray photons and the distance from the 
pulsar at which they are produced. For this purpose we modify 
the numerical code developed for the interaction and diffusion 
of electrons (Bednarek & Sitarek 2007). This code allows us not 
only to calculate the spectrum of y-rays produced by electrons 
but also determine their production sites around the pulsar, i.e. 
allowing us to study the morphology of the y-ray source. Since 
the electrons are immersed in a relatively strong magnetic field, 
especially close to the pulsar wind shock, we also include in the 
simulations their synchrotron energy losses during the diffusion 
process. We calculate the X-ray spectra produced by these elec- 
trons in the synchrotron process. In order to obtain reasonable 
precision of the IC y-ray spectra, we simulate the propagation 
of 1.5 X 10"* electrons per decade of the spectrum. The spectra 
are obtained within different regions around the pulsar defined 
by the radius /?Neb- 

We investigate the dependence of the X-ray and y-ray spec- 
tra on different parameters which determine the acceleration of 
the electrons (i.e. the magnetization parameter of the pulsar wind 
<T, the spectral index of the electrons' spectrum a\ the radius of 
the pulsar wind shock Rsh', and the minimum value of the mag- 
netic field in the nebula Bmin)- As shown in Fig. 2, the TeV y-ray 



spectra produced by the electrons in the IC process only weakly 
depend on the range of the considered parameters. On the other 
hand, the synchrotron X-ray emission depends on these param- 
eters much stronger (intensity, shape, energy range). The strong 
dependence of the synchrotron emission is due to the strong de- 
pendence of the magnetic field in the vicinity of the pulsar on the 
assumed parameters of the model. On the other hand, relatively 
weak dependence of the IC emission is due to the homogene- 
ity of the background radiation field (MBR and infrared galactic 
background) which is up-scattered by the relativistic electrons. 
We conclude that the TeV y-ray fluxes expected in this model de- 
pend rather weakly on the details of the electron spectrum (in the 
considered range of parameters). However, their intensity is ob- 
viously determined by the energy conversion efficiency from the 
pulsar to the relativistic particles. In contrast, the spectra of the 
synchrotron radiation in the X-ray range much stronger depend 
on the spectrum of the electrons and the propagation model. 

We also investigate the y-ray production in different volume 
around the Black Widow binary system B 1957+20. The IC y-ray 
and the synchrotron X-ray spectra are calculated assuming that 
this emission is produced within the region with the radius equal 
to 1.5 pc, 2.5 pc, 5 pc, 10 pc, and 15 pc (see Fig. 3). These dimen- 
sions correspond roughly to the angular size of the y-ray source 
on the sky equal to 2, 3.4, 7, 14, and 20 arc min for the distance 
of the source equal to 2.5 kpc. The electrons expand into such 
a region due to their diffusion in the nebula. Moreover, the TeV 
y-ray source is also expected to be shifted from the present loca- 
tion of the Black Widow binary due to its motion and/or limitted 
in specific directions by the diffusion of the electrons confined 
by the presence of the bow shock. In the case of a source with 
the radius above ~ 5 pc, the TeV y-ray source should appear ex- 
tended for the telescope array such as MAGIC. Our calculations 
show that most of the TeV y-ray emission (i.e. within a factor 
of two) is already produced within a region with the radius of 5 
pc. The shapes of the spectra, produced in specific parts of the 
y-ray source, are quite similar since the background radiation 
field (MBR and infrared), scattered by the relativistic electrons, 
fills this region homogeneously. Moreover the cooling process of 
the electrons is not very efficient. The electrons do not usually 
interact frequently but in a specific interaction lose significant 
amount of their energy when producing TeV y-rays. Due to the 
inefficient cooling, the parts of the spectra at low energies (in the 
GeV range), produced in the Thomson regime, are very similar. 
On the other hand, the synchrotron X-ray emission does not de- 
pend on the considered radius of the source at energies above a 
few keV. This can be understood since the hard synchrotron ra- 
diation is mainly produced close to the pulsar wind shock within 
the region with the extend of ~ 2 pc. There is however an im- 
portant contribution from the outer nebula to the part of the syn- 
chrotron spectrum at lower energies (below a few keV) since 
these electrons can still produce keV photons in the assumed 
minimum magnetic field. 

5. Comparison with observations of B1 957+20 

Finally, we compare the example calculations performed in 
terms of this modelling with the available observations of the 
Black Widow binary system B 1957 +20. The X-ray emission, ex- 
tending along the direction of the motion of the binary, has been 
detected by Chandra (Stappers et al. 2003, Huang et al. 2012). 
The X-ray synchrotron emission expected in our model has 
to be consistent with this observed spectral features. Recently, 
the pulsed GeV y-ray emission has been also reported from 
B1957+20 (Guillemot et al. 2012). The IC y-ray emission, pro- 
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Fig. 2. Gamma-ray (IC) and X-ray (synchrotron) spectra (Spectral Energy Distribution - SED) produced in the nebula around the 
Black Widow binary system containing the millisecond pulsar B1957-H20 for different model parameters. The spectra are produced 
by relativistic electrons which scatter the MBR and the infrared photons from the galactic disk. The maximum energies of the 
electrons are given by Eq. |4]and the minimum energies are equal \a - 1> TeV. (a) Dependence of SED on the magnetization 
parameter cr = 0. 1 (dashed), 0.01 (dotted), and 0.001 (solid) for the pulsar wind shock radius /?sh = 10'^ cm, the power law spectrum 
of the electrons with spectral index a - 2.5 and the minimum magnetic field strength Bmin - 0.5;uG. (b) Dependence of SED on 
the radius of the pulsar wind shock /?sh = 10'^' cm (dotted), lO'^ cm (solid), and 10'^ cm (dashed), for cr = 0.01, a - 2.5, and 
fimin - 0.5/iG. (c) Dependence of SED on the spectral index of the electrons a - 2.1 (dashed), 2.5 (solid), and 3 (dotted) for 
Rsh - 10'^ cm, cr - 0.01 and = 0.5/iG. (d) Dependence of SED on the minimum value of the magnetic field = 0.5/iG 
(dotted), 1/iG (solid), and 2yuG (dashed) for/?sh = lO'* cm, cr - 0.01, and a - 2.5. It is assumed that the relativistic electrons take 
10% of the rotational energy lost by the pulsar. The 100 hrs differential sensitivity of the MAGIC stereo system (thin dotted, Aleksic 
et al. 2012) and the 100 hrs CTA sensitivity (Actis et al. 201 1) are also marked. 



duced in the nebula by relativistic electrons, has to be below this 
pulsed emission. There are not available any positive detections 
or the upper limits on the TeV y-ray emission from this source. 

We have chosen intermediate parameters of the nebula from 
the range considered in Sect. 4. The IC and the synchrotron spec- 
tra are confronted with the available observations of the binary 
system containing B 1957+20 in Fig. 4. We have got good con- 
sistency with the level and shape of the X-ray spectrum from the 
nebula. Note that the X-ray observations put strong constraints 
on the parameters of the considered model. The emission extend- 
ing up to ~10 keV requires the presence of electrons with ener- 
gies at least ~ 4.7 x lO^B'^^J^ GeV (see Eq.©. On the other hand, 
the observed X-ray flux constrains the number of the relativistic 
electrons (which we fix on 10% of the pulsar energy loss rate) 
and the synchrotron energy loss rate which depends on oc B^EI. 
Therefore, it is not so easy to model the observed X-ray spectrum 



correctly since the change of the parameters have strong effect 
on the energies and intensity of the emitted synchrotron radiation 
(see calculations in Fig. 2). We conclude that observed X-ray 
extended emission put strong constraints on the parameters of 
the considered model. Having obtained consistency with the ob- 
served synchrotron spectrum, we calculate the IC y-ray spectrum 
for the same parameters (see the caption of Fig. 4). These spectra 
are confronted then with the sensitivities of the Cherenkov tele- 
scopes. We show the level of the IC emission expected from the 
region with the radius of ~5 pc, which is shifted from the present 
location of the binary system by about the same distance in the 
direction opposite to the movement of the binary due to the mo- 
tion of the binary system. Therefore, we conclude that the TeV 
y-ray source should be extended for the Cherenkov telescopes. 
The IC y-ray spectrum is clearly above the 100 hr sensitivity of 
the future Cherenkov telescope Array (CTA). It is also on the 
100 hr sensitivity limit of the MAGIC Cherenkov telescopes. 
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Fig. 3. SED of the y-ray spectrum from the IC process and the 
synchrotron X-ray emission produced by relativistic electrons in 
the nebula around the Black Widow binary system B 1957+20, 
integrated over a region with the radius: /?Neb = 1-5 (dot-dot- 
dashed curve) pc, 2.5 pc (dotted), 5 pc (dashed), 10 pc (solid), 
and 15 pc (dot-dashed). The other parameters of the model are: 
(J - 0.01, the pulsar wind shock radius 7?sh = 10'^ cm, the power 
law spectrum of the electrons with the spectral index a \ be- 
tween = 3 TeV and ^max (given by Eq.|4]i, and the minimum 
magnetic field strength B^in = 0.5;uG. The y-ray spectra are pro- 
duced by relativistic electrons which scatter the MBR and the 
infrared photons from the galactic disk. 



We conclude that even with the present Cherenkov telescopes 
(MAGIC, VERITAS) the bow shock nebula around the Black 
Widow millisecond pulsar B 1957 +20 might be detected. Note 
however that as the source is expected to be extended, the sensi- 
tivity of the present Cherenkov telescopes might become worse 
than the point source sensitivity shown in Figs. 2-4. 
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Fig. 4. The observations of the Black Widow binary containing 
millisecond pulsar B 1957+20: the X-ray tail emission detected 
by Chandra (Huang et al. 2012) and the pulsed, phase averaged 
y-ray emission discovered by Fermi (Guillemot et al. 2012) are 
compared with the calculations of the IC and the synchrotron 
emission by relativistic electrons in the nebula around this sys- 
tem. The y-ray spectra are produced by the electrons which scat- 
ter both the MBR and the infrared photons from the galactic 
disk. The electrons have a power law spectrum with the index 
of -2.1 between = 3 TeV and fimax - 160 TeV, given by 
Eq.|4] (thick dashed curves). The other parameters of the model 
are: the magnetization parameter cr = 0.01, the location of the 
shock /?sh = 10'^ cm, the minimum magnetic field strength 
fimin = 0.5/iG, the confinement region of the electrons has the 
radius of 7?Neb = 5 pc. The 100 hrs differential sensitivity of 
the MAGIC stereo system is marked by the thin dotted curve 
(Aleksic et al. 2012) and 100 hrs CTA sensitivity is marked by 
the thin dot-dashed curve (Actis et al. 201 1). 



6. Conclusion 

Assuming that the millisecond pulsars are able to accelerate 
electrons to relativistic energies in their vicinity, similarly as ob- 
served in the case of nebulae around classical pulsars, we calcu- 
late the synchrotron and the IC high energy emission from their 
nebulae. In fact, the existence of an extended synchrotron neb- 
ula has been recently confirmed in the Chandra observations in 
the case of the Black Widow binary system containing millisec- 
ond pulsar B1957+20. Therefore, as an example, we consider 
the bow shock nebula around this object. Note that in contrast 
to the nebulae around classical pulsars, the soft radiation field 
in the nebula around B 1957 +20 is not dominated by the syn- 
chrotron radiation but by the MBR and infrared radiation from 
the galactic disk. We have investigated the features of the X-ray 
and y-ray spectra for likely range of parameters which determine 
the nebula, assuming that the propagation of electrons is deter- 
mined by the diffusion process and/or the dynamical movement 
of the binary system. We conclude that the observed extended 
X-ray emission from the bow shock nebula can be explained by 
the synchrotron radiation of electrons provided that the energy 
conversion efficiency from the pulsar to the relativistic electrons 
is of the order of 10%. The TeV y-ray emission, produced by the 
same electrons in the IC scattering process, is expected to be de- 
tectable by the future CTA instrument. The predicted emission 



is also on the level of the 100 hr sensitivity limit of the MAGIC 
telescopes. However, since the nebula is expected to be extended, 
due to rather slow cooling process of electrons, the detectability 
of the TeV y-ray emission from the nebula around B 1957+20 
may be difficult. Note also that due to the motion of the binary 
system the TeV y-ray nebula should be shifted in respect to the 
direction towards the binary system by the distance comparable 
to the extend of the source (see also Cheng et al. 2006). 

Other bow shock nebulae around energetic pulsars should 
also emit synchrotron and IC high energy emission from their 
surrounding. However their detectability will strongly depend 
on the distance to the nebula. It can not be too large since 
the expected flux will be below detectability of the Cherenkov 
telescopes. But it should not be too close since the TeV y- 
ray nebula will have very large dimensions on the sky which 
again will make problematic its detectability with the Cherenkov 
telescopes. For example, the bow shock nebula around nearby 
Geminga pulsar (at the distance 169 pc) may not be detected 
by the present Cherenkov telescopes. Due to its small distance, 
the angular size of the TeV nebula expected in terms of dis- 
cussed above model, should be of the order of a few degrees, 
i.e more in accordance with the recent report on the presence of 
the extended ~20 TeV y-ray source with diameter (2.8±0.8)*', to- 
wards the Geminga pulsar by the MILAGRO observatory (Abdo 
et al. 2009). However, such nebulae might be detected by the 
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planned CTA which field of view can be as large as 8-9 degrees 
(Actis et al. 2011). 
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